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Abstract 
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OA ■ This letter derives the asymptotic symbol error rate (SER) and outage probability of multiple-input 



multiple-output (MIMO) maximum ratio combining (MRC) systems. We consider Rayleigh fading channels 



> 

o 

^ | \ with both transmit and receive spatial correlation. Our results are based on new asymptotic expressions 



which we derive for the p.d.f. and c.d.f. of the maximum eigenvalue of positive-definite quadratic forms 
in complex Gaussian matrices. We prove that spatial correlation does not affect the diversity order, but 
that it reduces the array gain and hence increases the SER in the high SNR regime. 
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I. Introduction 

Multiple-input multiple-output (MIMO) maximal ratio-combining (MRC) systems have recently received 
much attention due to their ability to mitigate the severe effects of fading through diversity [1]. The 
performance of MIMO-MRC has been previously investigated in various uncorrected and semi-correlated^ 
channel scenarios. The main performance measures evaluated have been symbol error rate (SER) and outage 
probability. Uncorrelated Rayleigh fading was considered in [2-7], and semi-correlated Rayleigh fading 
was considered in [8-10]. 

In this paper we present results for double-correlated channels. In particular, we derive new asymptotic 
expressions for the SER and outage probability for arbitrary numbers of antennas, which are simple 
functions of the system and channel parameters. These expressions are particularly useful in comparison 
to the only other known double-correlated MIMO-MRC results, presented in [11], which in any case had 
SER expressions limited to 2-antenna system^j . Moreover, the results we present here collapse trivially 
to semi-correlated scenarios, and provide simpler expressions than the previous results in [8-10]. 

Our results hinge on new first-order expansions which we derive for the probability density func- 
tion (p.d.f.) and cumulative distribution function (c.d.f.) of the maximum eigenvalue of positive-definite 
quadratic forms in complex Gaussian matrices. From these, we obtain a high SNR expression for the SER, 
and determine the diversity order and array gain. In so doing, we prove that MIMO-MRC achieves the 
maximum spatial diversity order, and we quantify the increase in SER in the high SNR regime due to 
spatial correlation. We also analyze the outage probability at outage levels of practical interest, and again 
quantify the effect of spatial correlation. 

II. MIMO-MRC System Model 
Consider an N r x N t MIMO system, where the N r x 1 received signal vector is 

r = vfHwi + n, (1) 



w 



where x is the transmitted symbol with E \x\ = 1, w is the beamforming (BF) vector with E 
1, n€ QN r y.i j s ac iditive noise vector ~ CNn t ,i (Ojv r xijIiV r )> an d 7 is the transmit SNR. Also, H is the 
N r x Nt channel matrix, assumed to be flat spatially-correlated Rayleigh fading. 

We assume that H can be decomposed according to the popular Kronecker correlation structure (as in 
[8, 10, 12]) as follows 

H = R5H W S2, (2) 



systems with two antennas at either the transmitter or the receiver. 



channels with transmit or receive correlation, but not both 
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where R > and S > are the receive and transmit spatial correlation matrices respectively, with unit 
diagonal entries, and H„, ~ CAfN r ,N t (ON r xN t ^N r ® I/v t ). 
The receiver employs the principle of MRC to give 

z = w+lTr = v^wWHwz + w^H^n. (3) 

Therefore, the SNR at the output of the combiner is easily derived as 

7 = ^wWHw. (4) 

The BF vector w is chosen to maximize this instantaneous output SNR, thereby minimizing the error 
probability. It is well known that the optimum BF vector w opt is the eigenvector corresponding to the 
maximum eigenvalue A max of H^H. In this case, the output SNR © becomes 

7 = 7 w lpt HtHw opt = T^max- (5) 

Clearly the performance of MIMO-MRC depends directly on the statistical properties of A max . For 
double-correlated Rayleigh fading channels, A max is statistically equivalent to the maximum eigenvalue of 
a positive-definite quadratic form in complex Gaussian random matrices Y ~ Q n ,m (In, (see [13] 

and [14] for more details), where m = max (N r , N t ), n = min (N r , N t ), and CI G C nxn and S G C mxm 
are the Hermitian positive-definite matrices 

A f R N r < N t A f S N r < N t 

n±{ v±{ (6) 

{ S N r > N t [ R N r > N t 

with eigenvalues uj\ < . . . < ui n and a\ < . . . < a m respectively. 

III. Maximum Eigenvalue Distribution of a Quadratic Form in Complex Gaussian 

Matrices 

We now present a theorem which gives new first-order expansions for the p.d.f. and c.d.f. of the maximum 
eigenvalue of a quadratic form in complex Gaussian matrices. This will allow us to derive the asymptotic 
distribution of the SNR in ©. 

Theorem 1: The following two expressions are first order expansions of the c.d.f. and p.d.f. respectively, 
of the maximum eigenvalue A max of Y ~ Q n ,m (In, 

F Am „(x) = ax mn + o{x mn ) (7) 



and 



/ Amax ( x ) = m nax mn - x + o (x" 1 "- 1 ) , 



(8) 



where 



a - 



T„ (n) 



det (n) m det(S) n r n (m + n) 
and T. (•) is the normalized complex multivariate gamma function defined as 



(9) 



r„(m) =JJr(m-»+l). (10) 

1=1 

Proof: The exact c.d.f. of the maximum eigenvalue of a positive-definite quadratic form in complex 
Gaussian matrices is given by [11, Theorem 1] 

(-l) n T n (n) det (ft)™" 1 det (S)™" 1 det (* (x)) 



^A max (x) 



(11) 



A n (f2)A m (£)(-xr (n - 1)/2 
where A m (•) is a Vandermonde determinant in the eigenvalues of the m-dimensional matrix argument, 
given by 



A m (£) = erf- 1 ] = 
Also, * (x) is an ra x m matrix with (i,j) th element 

where t = m — n, and 



n 



(12) 



P [m 



i ^ r 
i > t 



(13) 



2-1 



P(/;y) = l-e-^/A! 



(14) 



fc=0 



is the regularized lower incomplete gamma function. 

We seek a first-order expansion of (fTTT ) which will be immediate from a first order expansion of 
det (x)). Consider the Taylor expansion of det (x)) around the origin 



det (* (x)) = F (x) = F(<?) /°^ g + o (x Q ) . 



(15) 



g=0 



For a first order expansion we need to find the first non-zero coefficient in the sum. Using a well-known 
result for the qth derivative of a determinant, we have 

,1 /'/'' I*!/;},,/ 



{qu-,g m } 



ft! • • -9m! 



det 



(16) 



x=0 



where q\ + ■ — \- q m = q. Using ( fT3l ) and (fl4l) . we evaluate the derivatives in ( fT6l) as follows 



x=0 








i ^ t and (/i = 

i ^ t and gj > 

i > t and < qi < m 

i > t and qi > m 



(17) 



From (fTTT ). in order for the determinants in (fT6l ) to be nonzero, we require that for rows i ^ r we must 
have % = . We also require that for rows i > r, we must have qi > m, and these gfj's must all be 
different. Hence the smallest q for which these conditions are satisfied is given by 



q = m + (m + 1) + . . . + (m + (m — r — 1)) 
= mi! + n(n — l)/2. 



(18) 



For this smallest g value, we can now write (1161 ) as follows 



+ n) 



det 



x=0 



where the sum is over all permutations a = {a\, . . . , a n } of the numbers {0, . . . , n — 1}, and 



% < T 

m + Oj_ r i > r 



Using (fTTT ) in (fl9l we have 



T n (m + n) 



where 



'•j 



i 



m+Q!i_ T 



2 < T 



l> T 



(19) 



(20) 



(21) 



(22) 



We now focus on simplifying the determinant sum in (|2TI ). By removing the factors (l/o;j_ T ) ,n+cti ~ T from 
the determinants, and performing some row swaps, it can be shown that 

^ -J (_l)P-(«)(_l)r(r-l)/2 ^ M j 



Edet( S J=E(fl (;r- 

{q} {«} \i=T+l 

(_l)r(r-l)/2 det( (i 



det(n) m det(S) n ^ 11^. 



i=l 



(-ir (r - 1)/2 det((i-) 



det 



j'-i 



det(ft) m det(£) n 



(23) 
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where the last line followed from the definition of the determinant. Next, using ([T2l and the Vandermonde 
determinant identity [11, Eq. (56)] 



we can write ( f23b as follows 



0~j 0~ % 



lli=l°i 



Vdet (3 a ) = (- 1)-^- 1 )/ 2 (_ l)^-D/2 (_ Arn (^) A ra (fl 



det (S)" 1 ^" 1 det (ft)"^" 1 



_ / ^ntm+i) A m (S) A n (fl) 

1 J det(S) m+n - 1 det(fi) m+n - 1 

Substituting (1251) into (f2TT) and using (TT8T ). we derive the desired first-order expansion of det (x)) in 
(fl"5l) . given by 



(lr , ( vTr [ , ) ) - (-l) n (-l) n( " 1)/2 A m (S) (fl) +n(n _ 1)/2 mn+n(rl _ 1)/2 

r ra (m + n) det (2J) det (&2) 



The c.d.f. result (O now follows by substituting d26l ) into (fTTT ) and simplifying. The p.d.f. result ([8]) then 
follows trivially by taking the derivative of (0 w.r.t. x. □ 
For the special case of uncorrected fading, (O reduces to 

F (n) 

Fx max (x) = = "V , x mn + o (x mn ) . (27) 
r n (m + n) 

which agrees with a result derived previously in [7]. 

IV. Asymptotic SER and Outage Analysis of MIMO-MRC 
For many general modulation formats, the average SER of MIMO-MRC can be expressed as 

P s = £ 7 \oQ (\p%n 



(28) 



where Q (•) is the Gaussian Q-function, and a and b are modulation-specific constants [15]. In [16], an 
alternative expression of (|28T ) was provided as follows 

a\/b f' 00 e~ bu 

Ps = ^ e -^F x _{u)du. (29) 

We now analyze the SER performance in the high SNR regime in order to derive the diversity order 
and array gain of the system. Armed with Theorem [TJ we can directly invoke a general parameterized 
single-input single-output (SISO) SER result from [17], and perform some basic algebraic manipulations, 
to obtain a high SNR SER expression given by 

SER°° = (G a • 7)~ Gd + o {T Gd ) , (30) 
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where the diversity order is 



Gd = mn 



(31) 



and the array gain is 



G a = det (ft) 1/n det(£) 1/m 26 



aT n (n) 



) 



l/mn 



2T n (m + n) 



(2mn - 1)!! 



(32) 



with 



(2mn — 1)!! = 1 x 3 x . . . x (2mn — 1) . 



(33) 



We clearly see that MIMO-MRC achieves the full spatial diversity order of mn, regardless of the spatial 
correlation. Moreover, using Hadamard's inequality and the fact that the diagonal elements of ft and £ 
are unity, it is easily found that 



with equality in the upper limit only when the correlation matrices are identity matrices. Hence, from (|32~1) 
we see that the effect of the correlation is to reduce the array gain (with respect to uncorrected fading) 
by a factor of det (ft) 1/n det (£) 1/m , thereby increasing the SER in the high SNR regime. Note that for 
the special case n = 2, this result can be shown to reduce to an expression reported previously in [11]. 

We now consider the outage probability of MIMO-MRC systems in double-correlated Rayleigh channels. 
The outage probability is an important quality of service measure, defined as the probability that 7 drops 
below an acceptable SNR threshold 7^. It is obtained using ((5]) as follows 



In practice, we are usually interested in small outage probabilities (i.e. 0.01, 0.001, ...), which correspond 
to small values of 7th. To gain further intuition at these small outage probabilities, we use CO) in Theorem 
[7] to write the outage probability in (I35T ) as follows 



According to d34l) , this result shows explicitly that in the low outage regime the outage performance 
degrades due to the presence of spatial correlation. Moreover, the increase in outage probability (with 
respect to uncorrected fading) is quantified by the factor det (fi)~ m det (E)~ rt . 



sC det (ft) ^ 1 and < det (S) < 1 



(34) 




(35) 




(36) 
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V. Numerical Results 

For our numerical results we construct the correlation matrices using the exponential correlation model. 
The (i,j)th entries of Q and S are given by {O}^ ^ = p± and {S}^- = p\ with p\,p2 G [0,1), 
respectively. Note however, that all of the analytical results presented in this paper apply equally to any 
other correlation model that conforms to the general structure in (f2]). 

Fig. □ shows the SER of MIMO-MRC with 8PSK (a = 2, b = 0.146) modulation, for various antenna 
configurations. The 'Analytical' curves are generated via numerical integration of d29l ) using the analytical 
c.d.f. (TTTb - The 'Analytical (High SNR)' curves are based on (l30b . Clearly the diversity orders and array 
gains predicted by the high SNR analytical results are accurate. Also, as expected from the analysis in 
Section |IVJ we see that the SER increases monotonically with the level of correlation for both antenna 
configurations. 

Fig. |2] shows analytical and Monte-Carlo simulation outage probability curves for a MIMO-MRC system, 
comparing different correlation scenarios. The analytical results are based on (TTTb . As expected from our 
asymptotic analysis in Section JVJ we see that the correlation increases the outage probability for all outage 
levels of practical interest (i.e. in this case, for outage levels < 30%). It is also interesting to observe that 
the opposite occurs for high outage levels. 

VI. Conclusions 

We have examined the asymptotic performance of MIMO-MRC systems in double-correlated Rayleigh 
channels. Our results are based on new closed-form asymptotic expressions which we have derived for 
the marginal maximum eigenvalue distribution of positive-definite quadratic forms in complex Gaussian 
matrices. The new results prove that the presence of spatial correlation yields a net increase in SER in the 
high SNR regime, and also degrades the outage performance for outage levels of practical interest. 
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SNR(dB) 

Fig. 1. Symbol Error Rate of MIMO-MRC in various double-correlated Rayleigh channels with 8PSK 
modulation. 




Fig. 
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2. Outage probability of 3 x 3 MIMO-MRC in various double-correlated Rayleigh channels, and for 
OdB. 



